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Utilizing a full microscopic Bogoliubov-de Gennes (BdG) approach, we study the equilibrium charge and 
spin currents in ballistic SFSFS Josephson systems, where f is a uniformly magnetized ferromagnet and 
S is a conventional ,s-wave superconductor. From the spatially varying spin currents, we also calculate the 
associated equilibrium spin transfer torques. Through variations in the relative phase differences between the 
three S regions, and magnetization orientations of the ferromagnets, our study demonstrates tunability and 
controllability of the spin and charge supercurrents. The spin transfer torques are shown to reveal details of the 
proximity effects that play a crucial role in these types of hybrid systems. The proposed S FS FS nanostructure 
is discussed within the context of a superconducting magnetic torque transistor. 
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I. INTRODUCTION 

Proximity effects inherent to superconducting systems with 
inhomogeneous magnetic order presents a mechanism by 
which dissipationless current flow and the spin degree-of- 
freedom can both be effectively coupled and controlled.^® 
The important role that proximity effects play in the static 
and transport properties of ferromagnetic Josephson junctions 
with i-wave superconductors is now well established. In¬ 
deed, the proximity induced damped oscillatory supercon¬ 
ducting correlations within the ferromagnet region serves as 
a channel for interlayer coupling and spin switchin g^ P | 5 | 7 | 9H ii] 
Proximity induced triplet pairing correlations within the fer¬ 
romagnetic Junction also provides another avenue for spin 
transport.EEH Interest in Josephson junctions with ferromag¬ 
netic layers has grown due to their possibility as serving as 
elements in next generation superconducting computing and 
nonvolatile memories!^ where single flux quantum circuits 
containing multiple Josephson junction arrangements can im¬ 
prove switching speeAs.^^G^To determine whether Josephson 
structures can serve as viable cryogenic spintronic devices, it 
is crucial to understand the behavior of the spin currents that 
can flow in such systems. The spin current flowing into the 
ferromagnetic regions exerts a torque on the magnetization if 
the current polarization direction is noncollinear to the local 
magnetization in the ferromagnet. In other words, the spin 
angular momentum of the polarized current w ill be partially 
transferred to the magnetization in the F region.EEHlThis spin 
transfer torque (STT) serves as an important mechanism in 
spintronics devices.n^El The STT effect can cause magne¬ 
tization switching for sufficiently large currents without the 
need for an external held. This switching aspect provides a 
unique opportunity to create and improve fast-switching mag¬ 
netic random access memories.E2^ 

The recent experimental pursuits of spin-based memory 
technologies involving various arrangements of S FS Joseph¬ 
son junctions has rekindle d in terest in the realm of ferro¬ 
magnetic Josephson arrays.l3®tS] When a sequence of S FS 
junctions are placed in a series configuration, creating a 
SFSFS type junction shown in Fig. additional possi¬ 


bilities emerge for the control of the associated spin and 
charge supercurrents.1212^ For example, the triplet components 
of the supercurrent and total charge transport in diffusive 
S FS FS structures is closely linked to the relative magneti¬ 
zation orientations, which can directly alter the total charge 
current flow, causing it to reverse direction in the ferromag¬ 
netic layers .1^ The transport of triplet supercurrents through 
the middle S electrode can be utilized to manipulate the mag¬ 
netic moment of the F layers in SFSFS hybrids.!® The spin- 
polarized supercurrents in these types of systems may also 
be used to indu ce a STT acting on the magnetization of a 

ferromagnet.l2122J!lll 

At the interfaces between the F and S regions in a ballis¬ 
tic S FS FS Josephson junction, q uasip articles undergo An¬ 
dreev and conventional reflections.!®^ Besides the contribu¬ 
tion from the continuum states, the superposition and interfer¬ 
ence of the corresponding quasiparticle wavefunctions in the 
F regions result in subgap bound states that contribute to the 
total current flow between the S banks. By varying the width 
of the central S layer, ds , modifications to the Andreev bound 
state spectra can ensue, i.e., by simply decreasing d^, addi¬ 
tional overlap can occur between the subgap bound states in 
the adjacent F regions .!®02 For S FS FS type structures, if the 
central S layer is sufficiently thin, i.e., ds < the proximity 
effects within the interacting F regions result in the Cooper 
pair amplitude and local density of states in each F region be¬ 
ing mutually altered. Therefore, magnetization rotation in a 
single F layer can strongly influence the thermodynamic and 
quantum transport properties throughout the rest of the sys¬ 
tem. The coupling between the different regions can then re¬ 
sult in the system residing in abound state corresponding to 
a phase difference of Aip - The appearance of superhar¬ 
monic Josephson currents (with second and higher harmon¬ 
ics: sin2(^, sin 3(f, ...) were theoretically predicted to appear 
in nonequilibrium and point contact Josephson junctionJ®. 
Shortly thereafter, the higher harmonic supercurrents were ex¬ 
perimentally observed in nonequilibrium situation^®. Also, it 
was shown theoretically for a uniform SFS junction, that the 
higher harmonics can be revealed at the O-tt transition point, 
where the first harmonic is highly suppressed due to the super¬ 
current flow reversing direction at that point.!® Therefore, the 
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nonvanishing supercurrent at the O-tt transition point observed 
experimentally in R ef.l45]wa s soon attributed to the presence 
of higher harmonicJ3^®^. Subsequent works with ferro¬ 
magnetic Josephson junctions demonstrated that the higher 
harmonics can naturally arise when varying the location of 
domain wallJSI, and in ballistic double magnetic SFFS junc¬ 
tions, provided that the thickness of the magnetic layers are 
unequaP. Recently, evidence of higher harmonics has been 
experimentally observed in Josephson junctions with spin de¬ 
pendent tunneling barriers.® 

The focus of this paper is to theoretically investigate prox¬ 
imity effects leading to modified superconducting correlations 
and controlled charge and spin transport in SFSFS ballis¬ 
tic junctions. We will address a variety of relative magne¬ 
tization orientations, and prescribed superconducting macro¬ 
scopic phase differences between the S terminals. Utilizing a 
microscopic Bogoliubov-de Gennes (BdG) approach, we de¬ 
rive the appropriate expressions for the charge and spin cur¬ 
rents and the corresponding equilibrium spin transfer torques. 
The numerical solutions to the BdG equations are employed 
to study the current phase relations, revealing the emergence 
of additional harmonics that depend on the tunable magneti¬ 
zation profile and other system parameters. We demonstrate 
that our proposed SFSFS systems can be considered as a 
superconducting magnetic torque transistor, where the flow 
of spin and charge currents can be tuned by the macroscopic 
phases of the superconducting leads. This, in turn, dictates the 
torques acting on the exchange fields of the F layers. Remark¬ 
ably, the superconducting phases (in addition to other system 
parameters) can effectively switch the torques acting on the 
magnetizations of the F layers ‘on’ or ‘off’. The directions 
of the torques and charge currents are shown to not be re¬ 
lated by simple functions of the phase differences or exchange 
fields, similar to what was observed in simpler ferromagnetic 
Josephson junctions involving equilibrium torques.l^ We also 
find that when the angle describing the relative in-plane mag¬ 
netic exchange field orientations is varied, the torque tending 
to align the two F magnetizations is usually largest for rela¬ 
tive magnetization angles other than the expected orthogonal 
configurations. Moreover, we have found that these sequential 
nanodevices allow for detecting pure second harmonics in the 
current phase relations, depending on the system parameters, 
including the relative magnetization orientations. We present 
a study of the crossover between the first and second harmonic 
in the current phase relations and consider experimentally fea¬ 
sible situations to observe them. This crossover is discussed 
in the context of the appearance of equal spin triplet correla¬ 
tions with m = +1 spin projections along the spin quantization 
axis. 


The paper is organized as follows. In Sec.|I^ we outline the 
theoretical approach used and derivation of various physical 
quantities investigated, including the supercurrent, magneti¬ 
zation, spin current, and associated torques. We present our 
results in Sec. [Ill] This section is divided into two-subsections: 
Subsect. Ill A presents the current phase relations and the sec¬ 
ond harmonic supercurrents that can be generated by calibrat¬ 
ing the system parameters. Subsect. |IIIB discusses the as¬ 
sociated spin currents and equilibrium spin transfer torques. 



FIG. 1. (Color online) Schematic of the SFSFS ballistic hetero¬ 
junction considered in this paper. The system is infinite in the 
yz plane, and thus the x axis is normal to the plane of the inter¬ 
faces. The middle S electrode has a thickness dg in the x direc¬ 
tion, and is sandwiched between two F layers with unequal thick¬ 
nesses dfi and df 2 - The F layers are uniformly magnetized and 
their exchange fields are denoted by h, (for i = 1,2). To sim¬ 
plify notation, we have defined the magnetization directions via: 
hi = l/r/Kcosa,, sin a,- sin/?,, sin or, cos/l,). The S electrodes can take 
arbitrary phases, described by <pi, ipM, and ipR for the left, middle and 
right S electrodes, respectively. 


Finally, we give concluding remarks in Sec. IV 


II. THEORETICAL METHOD 


We begin our methodologyby introducing the Bogoliubov- 
de Gennes (BdG) formalism.ESThe BdG approach is a conve¬ 
nient microscopic quantum mechanical technique that allows 
a complete investigation into the fundamental characteristics 
of the superconductivity of ballistic superconducting hetero¬ 
junctions. The microscopic BdG formalism can easily accom¬ 
modate a broad range of magnetic exchange field strengths 
and profiles, including the half-metallic limit where the mag¬ 
nitude of the exchange field and the Fermi energy, bf, are the 
samd^. A schematic of the multilayer configuration that we 
study is depicted in Fig. For this quasi one-dimensional 
system, physical quantities are invariant with respect to the 
yz plane, while the x-direction captures the essential phys¬ 
ical characteristics of the system. The corresponding spin- 
dependent BdG equations are thus expressed as. 


FIq - —hx + ihy 0 

-hx — ihy Ffo + A(x) 

0 A*(x) -CKo - K) 

A*(x) 0 —hx + ihy 


A(x) 

0 

—hx - ihy 
-(■^ 0 +/*.). 




= en^n(x), 


( 1 ) 


where T'„(x) = (m„|(x), m„j,(x), v„t(x), v„j.(x))’', and u„a- and 
v„a- are the quasiparticle and quasihole amplitudes. The pair 
potential A(x), which effectively scatters electrons into holes 
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(and vice versa) is nonzero only in the superconducting elec¬ 
trode regions. We furthermore assume that A(x) is piece- 
wise constant in the S regions, with each S region possess¬ 
ing the same magnitude but possibly different phase. Thus, 
within the external S electrodes, A(x) takes the form 
in the left, Aoe"^" in the middle, and Aoe"^® in the right S 
electrode. The combinations of phase differences involving 
ifL, (fM, and ifR results in additional possibilities for super¬ 
current flow compared to conventional Josephson junctions 
comprised of two superconducting banks. The single parti¬ 
cle Hamiltonian '7fo(x) is defined as, 

1 

'f^oix)^- — -^+e±-Sf + U(x), ( 2 ) 

Im ox^ 

where e_i - ^(ky + k^) is the quasiparticle energy for motion 
in the invariant yz plane (see Fig.[T]l, and the spin-independent 
scattering potential is denoted by U(x). We represent the mag¬ 
netism of the F layers by a Stoner effective exchange energy 
h{x) which will in general have components in all {x,y,z) 
directions. Additional technical details on solving the BdG 
equations for this type of quasi one-dimensional setup is given 
in the Appendix. 

Various other types of “inverse” proximity effect^HEl 

can 

also occur in the vicinity of the F/S contacts, whereby ferro¬ 
magnetic order propagates from one F layer to the other (cre¬ 
ating a mutual torque) via the central S layer. Therefore it is of 
interest to determine not only the spatial profile of the magne¬ 
tization m(x) within the F regions, but also within the central 
S layer, where the induced magnetization can also screeiJSI 
the magnetization of the adjacent ferromagnet. The complete 
spatial profiles of the magnetization are determined using the 
expressions!^ 


in terms of the quasiparticle amplitudes !®3 


Jxix) = — ^ Im j /„ + Uni \ 


dx dx 

dv\ dv* 


r ov* civ* A 

+ ]J. 


Similarly, following the approach outlined in the Appendix, 
we can also write the spin current So- with spin cr flowing 
along the x direction in terms of the quasiparticle amplitudes: 
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-d -.4^1}. (9) 

where the sums for the currents above are in principle taken 
over all quasiparticle states. 


III. RESULTS AND DISCUSSIONS 


nixix) = -fiB '^[{uni{x)Uni{x) + M*^(x)M„t(x))/„ 
n 

Iflyix) = -i^lB^[{u„i(x)ul^(x) - M„i(x)M*^(x))/„ 
n 

+ [v„|(x)v*^(x) - v„|(x)v*^(x)](l - /„)}. (4) 

m^(x) = -jUij^([|M„|(x)|^ - \Uni(x)f-]f„ 
n 

+ [|v„TWl'-|v„i(^)l'](l-/4 (5) 

where is the Bohr magneton, and f„ the Fermi function. 

The Josephson effect leads to many possibilities for charge 
supercurrent transport in SFSFS junctions. This experimen¬ 
tally accessible phenomenon is now well understood, with the 
primary driving mechanism being the difference between the 
mac roscopi c phases of two S banks, separated by a weak 
link.l^l— When there are three coupled S banks, the sit¬ 
uation becomes more complicated, and the dissipationless 
charge current depends on various combinations of the phase 
differences in addition to the other geometric and material 
properties of the system. When computing the supercurrent 
flowing in the x direction, we express the Josephson current 


We focus here on the low temperature regime, with T jTc - 
0 . 001 , where Tc is the critical temperature of the correspond¬ 
ing bulk S material. For simplicity we set cpi - 0, ifM - <,p/2, 
and (fB - (f, for the phases of the left, middle, and right S ter¬ 
minals, respectively. Thus, a phase difference of (^/2 is main¬ 
tained across each S electrode. The spatial variables are nor¬ 
malized in terms of the Fermi wavevector, including the BCS 
zero-temperature coherence length, set to kf^o = 100 , and 
the dimensionless position X, written as A = kfX. For each 
physical quantity studied, a broad range of central S widths 
will be considered. We assume that the ferromagnets are 
similar materials with identical exchange held strengths, i.e., 
IJiil = \h 2 \ = h, set to the representative value of /r/e/r = 0 . 1 . 
To create favorable conditions for equal-spin triplet genera- 
tion,l^ the Fi and F 2 regions have highly asymmetric widths, 
with dfi - O.l^o, and df 2 - 3. 8 ^ 0 , so that dfi dF 2 - 

A. Josephson Charge Supercurrent 

We begin with a discussion of the supercurrent charge trans¬ 
port by solving the microscopic BdG equations (Eq. Q) over 
a broad range of energies and then summing the correspond¬ 
ing quasiparticle amplitudes and energies according to the ex¬ 
pression given by Eq. The charge current also can be 
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FIG. 2. Normalized Josephson current versus the macroscopic phase 
difference of the two outermost electrodes, A^, in a 5 FS FS struc¬ 
ture with asymmetric ferromagnet widths of = 0.1 and 

dfil^o = 3.8 (see Fig. 0 The width of the central S electrode, dg , 
varies as shown in the legend. The relative exchange fields between 
the two magnets is orthogonal with ai = a 2 = 90°, Pi = 90° (along 
y), and p 2 = 0° (along z). The phase of the middle 5 electrode takes 
the value, <Pm = i^/2. For comparison, in (a) we show the results for 
a simpler 5 FS junction having a phase difference Aip/l and F width 
dpi = 0 . 1 ^ 0 . 


obtained by minimizing the free energy with respect to the 
appropriate superconducting phase differences.E^Our micro¬ 
scopic method fully accounts for bound states that may be 
generated from quasiparticle trajectories with large momenta 
corresponding to in-plane energies comparable to ef (shown 
to be important for SNS junction^^. The charge supercur¬ 
rent is normalized by Jq - nevp, where vp is the Fermi ve¬ 
locity, e the electron charge, and n is the number density. We 
focus our attention on supercurrents flowing through the fer- 
romagnets, recalling that each of the three S regions act as ef¬ 
fective sources or sinks in the current. The pair potential A(.r) 
of course vanishes in the intrinsically nonsuperconducting F 
regions. 

To begin, in Fig. the supercurrent is shown as a function 
of the phase difference, A^, for a wide range of central S elec¬ 
trode widths, ds ■ The central S electrode acts as an external 
current source, and hence the spatial behavior of the current is 
piecewise constant in each F region. Thus, each panel corre¬ 
sponds to the current in a particular ferromagnet (as labeled). 
The relative magnetic exchange fields are orthogonal, with hi 
directed along y and h ,2 along z (see Fig. [^. Two limiting 
cases are shown: In the first case, the width of the central S 
layer is zero {dg = 0 ), and in the second case, a large cen¬ 
tral layer (dg - 5 ^q) is considered. When there is no middle 
S layer, the current phase relation (CPR) is 7 r-periodic, with 
behavior consistent with a ballistic S FFS asymmetric double 
magnetic structur^. For ds = 5^o, the large S width effec¬ 
tively decouples the two ferromagnets, creating two isolated 
SFS junctions with phase differences ip 12. Thus, in this case 
one junction consists of a thin uniform Fi region sandwiched 
by two superconductors, and its CPR reflects the overall be¬ 
havior and direction reversal that is emected in narrow fer¬ 
romagnetic SFS Josephson junctions.™ The other decoupled 
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FIG. 3. Normalized supercurrent flowing through the ferromagnetic 
regions vs Ap. The central S layer has width dg = ^o- We con¬ 
sider several magnetization orientations, Pi (see legend), inside the 
Fi layer of width dpi = O.lfo- The magnetization direction of the 
larger F 2 layer {dp 2 = 3. 8 ^ 0 ) is strictly along z, corresponding to 
P 2 = 0 °. 


SFS junction containing F2 (of width dF2 - 3. 8 ^ 0 ) has a sub¬ 
stantially diminished current due to its much greater width. 
Note that there are negligible triplet correlations present when 
ds IS large due to an effectively uniform magnetization in each 
F layer. For intermediate S layers, the CPR evolves from its 
form in one of these limiting cases, to a richer more complex 
one due to the emergence of additional harmonics. This is 
due in part to the greater amount of triplet correlations that 
are present when the F layers possess orthogonal magnetiza¬ 
tion configurations. The appearance of additional harmonics 
in the current phase relation has been discussed in the diffu¬ 
sive and clean reg imes for simpler ferromagnetic Josephson 
junction structures.^SElllIl Note that in Fig. |^a), the phase 
corresponding to the hrst peak in the current phase relation, 
denoted by the critical phase, p*, has p* ^ 60° for ds - 0 , and 
then gets smaller for ds ^ ^0 ? before increasing nearly lin¬ 
early with ds ■ This is contrast to what is observed in Fig.j^b), 
where p* increases monotonically with ds ■ 

Figure 1^ shows the CPRs at various magnetization direc¬ 
tions Pi. The magnetization in F 2 is fixed along the z direc¬ 
tion. The width of the central S layer is now set at ds — 
and as stated earlier, its phase has the value pjl to ensure that 
adjacent superconductors maintain the same phase difference. 
Examining panel (a), we see that when the relative magneti¬ 
zations are parallel (fii - 0°) or antiparallel {J3i - 180°), the 
current exhibits a nearly sinusoidal CPR. For intermediate Pi 
leading to noncollinear magnetizations, higher order harmon¬ 
ics appear in the Josephson current. Figure [^b) reveals that 
in the wider F2 region, collinear orientations result in regu¬ 
lar sawtooth-like patterns in the charge current as Ap varies. 
The current in the larger magnet flows in opposite directions 
depending on whether the relative magnetizations are parallel 
or antiparallel, in contrast to the narrow Fi segment reported 
in (a). Similarly to what is observed in the narrow Fi region, 
we also find more complicated higher order harmonics in de¬ 
scribing the current for misaligned relative magnetizations. 
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FIG. 4. (a) and (b): Normalized Josephson current vs the relative 
in-plane magnetization angle, P\ (see Fig.[TJ. The exchange fields 
in the two magnets are parallel when Pi = 0, and antiparallel when 
Pi = 180°. The magnitude of the exchange field is set to h = O.ls/r. 
The outer S electrodes have A<^ = 45°. Five different S widths are 
considered, as depicted in the legend. Panels (c) and (d) show the 
average (over the central S region) of the magnitudes of the triplet 
correlations vs jSj. 


In the broader context of layered FjS structures, includ¬ 
ing ferromagnetic Josephson junctions and spin valves, mis¬ 
alignment of adjacent F layer magnetizations will typically 
generate equal-s pin p airing that is greatest in the orthogo¬ 
nal conhguratiorP^l®] Jq investigate the supercurrent trans¬ 
port properties when proximity-induced triplet pair correla¬ 
tions are present in S FS FS type structures, it is instructive to 
investigate the sensitivity of to relative magnetizations ori¬ 
entation. Therefore in Fig.j^ the Josephson current and triplet 
correlations are shown in the two ferromagnets as a function 
of exchange field orientation, Pi. As yfij sweeps between the 
parallel (ySi = 0°) and antiparallel (ySi = 180°) states, adjacent 
S electrodes are, as before, maintained at constant phase dif¬ 
ference ipjl. Multiple middle S terminal thicknesses are con¬ 
sidered (see legend), with the dg -0 curve shown for compar¬ 
ison purposes. As shown in Fig.Qa) and (b), when magnetic 
coupling is significant (for ds < 2^o) the supercurrent is non¬ 
monotonic and can be highly sensitive to the relative direction 
of the magnetic moments in the F layers. When the mag¬ 
nets have collinear magnetizations, the Josephson current in 


Fi is often weaker for the parallel configuration compared to 
the antiparallel configuration, except when there is no central 
superconductor, or when it is very wide. The orthogonal state 
(jSi = 90°) however results in the maximal current flow. Incre¬ 
ments in the central electrode thickness can drastically mod¬ 
ify the supercurrent signature. Eventually for large enough 
ds, the magnetic coupling is diminished, and variations in Pi 
can no longer affect the current flow. The corresponding de¬ 
coupled S FS junctions then have uniform current flow, that is 
larger in the narrow Fi region (panel (a)), and in Fi (panel(b)), 
becomes negligible due to the larger width. 

The emergence of additional harmonics in the current phase 
relations is often correlated with the generation of triplet cor¬ 
relations that are odd in timJSlor frequency. As we saw previ¬ 
ously in Fig. varying the phase difference A(f, revealed the 
emergence of additional harmonics as the relative exchange 
fields went from the parallel {P 2 = 0°) to orthogonal {P 2 = 
90°) magnetic configuration. To further explore the evolu¬ 
tion of triplet pairing correlations with magnetic orientations, 
in (c) and (d) the spatially averaged triplet amplitudes |/o,avgl 
(with spin projection m = 0), and |/i,avgl (with spin projection 
m - +1) are shown as functions of Pi. These quantities are 
calculated using the expressions ^3/o(v, t) - l/22„(/J'*'(v) - 
and/i(x,f) = 1/2 where 

we define („(t) = cos(e„f) - isin(e„f)tanh(e„/(27’)), and 
= Mno-(v)v*^,(v). The summations are in principle over 
all states. A representative value of f = 6 is used for the scaled 
relative time, where I = cjot. The quantization axis in the re¬ 
gions of interest is aligned along the z-direction, however it 
is straightforward to align it along a different axis that may 
coincide with the local magnetization direction.El Comparing 
Figs.j^c) and (d), it is evident that the behavior of the triplet 
amplitudes as a function of Pi is anticorrelated, with the aver¬ 
age 1 / 0 1 smallest when |/i| peaks at /Ji = 90°. Increasing the 
S width is shown to reduce the /o amplitudes gradually, how¬ 
ever the equal-spin component f\ drops much more abruptly 
to negligible values once ds exceeds ^o- Although not shown, 
the singlet correlations within S (for all S widths) were found 
to not exhibit significant sensitivity to changes in P\. This is 
clearly in sharp contrast to what is observed for both triplet 
components. Having discussed now some salient features of 
the charge currents, we now turn our attention to spin trans¬ 
port and the corresponding equilibrium spin transfer torques 
within the junction region. 


B. Spin currents and spin transfer torques 

Spin-polarized transport quantities are of paramount impor¬ 
tance when studying S FS FS type junctions as potential com¬ 
ponents in spintronics devices. The spin current S' is a local 
quantity responsible for the change in magnetizations due to 
the flowing of spin-polarized currents. The main contribu¬ 
tor to the equilibrium spin current and corresponding spin- 
transfer torque t is the spin-resolved Andreev bound stateJSI 
which play the main role in torque sensitivity to variations in 
A(f and Pi. Thus, the STT can be a useful probe of the spin de¬ 
gree of freedom in S/F proximity elements. The current that 
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is generated from the macroscopic ph ase diffe rences in the S 
electrodes can become spin-polarizecP®^!^ when entering 
one of the fetTomagnet regions. A portion of this spin current 
can then interact with the other ferromagnet and be absorbed 
by the local magnetization due to the spin-exchange interac- 
tions.lSl Since we are considering ferromagnets with in-plane 
magnetic exchange fields, the only spin current that can flow 
is the out-of-plane component Sx- This is consistent with the 
fact that only Tx can exist in equilibrium when spin currents 
do not enter or leave the superconducting electrodes.® 

As shown in the Appendix, the method used here to de¬ 
termine T involves simply calculating the magnetic moment 
throughout the entire system and then using. 


I^B 


-m X h. 


( 10 ) 


where the magnetization components are given in Eqs. 0 - 
0 - Equivalently, in the steady state, one can use the conti¬ 
nuity equation for the spin current [Eq. (B31] to determine the 
torque transferred by simply evaluating the derivative of the 
spin current as a function of position: 


T — 


dx ' 


( 11 ) 


The net flux of spin current AS x through a certain region 
bound by points xi and X 2 is therefore: 


rx2 

AS x^ Sx{X 2 ) - Sxixi) ^ I dxTx^Tx,iot- (12) 

Jx, 


In other words, the change in spin current at the interface 
boundaries (x - xi and x - X 2 ) is equivalent to the net 
torque acting within those boundaries. Either approach, us¬ 
ing Eq. ([Tg, or Eq. ( [TT] i, is sufficient to calculate r, as they 
both yield precisely the same result. In the results that fol¬ 
low, we calculate the torques using Eq. ( [T0| ), thus avoiding the 
numerical derivatives that arise when using Eq. 0 . 

We first present in Eig. |^a) the .^-component of the local 
spin curi'ent, Sx, normalized by So = -hbNfSfI kp, where 
Nf is the density of states at the Eermi energy. We numer¬ 
ically calculate 5;^ by summing the quasiparticle amplitudes 
and energies using Eq. (0. Several different phase differences 
are studied as shown in the legend, and the exchange interac¬ 
tions are orthogonal: h\ — {Q,hyi,Q), and h .2 = (0,0, h^ 2 )- 
The central S layers is one ^0 wide. The spin cutTent reveals 
precise spatial behavior of the junction interlayer magnetic 
coupling, and from Eq. (Ill, one can deduce the correspond¬ 
ing local behavior of t^. In F 2 , the oscillating spin currents 
each have a phase and magnitude that can change, depending 
on Aip. Once the spin current enters the S region (bound by 
the dashed vertical lines), it immediately becomes conserved, 
whereby there is no transfer of spin angular momentum. Since 
dS XI dx = 0, we have = 0 in that region. Within Fi, the nar¬ 
row width limits the extent at which S x can vary, and conse¬ 
quently it undergoes a nearly monotonic decline, before van¬ 
ishing within the superconductor. Since 5;^ = 0 in the outer S 
electrodes, the difference AS x over either ferromagnet is de¬ 
termined by the value of 5 ^ within the central S. Thus, despite 



A(p A(p 


FIG. 5. (a) Normalized spin current, Sx, as a function of normalized 
position X. Several phase differences A(p are considered. The in¬ 
terfaces separating each region are denoted by vertical dashed lines. 
The spin current is conserved in the central S region, where there is 
no magnetic exchange interaction, (b) Normalized total equilibrium 
torque, Tx,tot, as a function of the superconducting phase difference 
A(fi between the outermost S electrodes. A wide range of dg /fo ra¬ 
tios are considered (see legend). 


drastically different local behavior of 5 x in each F region, the 
net flux of spin current through either Fi or F 2 differs only in 
sign. One can then see that by examining the value of the con¬ 
served 5x in the central S region, the flux ASx through e.g., 
El is largest when Aip - 180°, and smallest when Atp = 120°. 
This observation is consistent with (b), where the total torque 
(normalized by tq = -iibNfEf) is shown as a function of 
Aip. We calculate Tx,tot over the F\ region using Eq. (12i, 
although the result for F 2 is trivially obtained, since within 
each of the three i-wave S electrodes, there can be no flux of 
spin current. This requires Tx,tot in Fi to be the exact opposite 
in E 2 . As seen in (b), the net torque can be quite sensitive 
to the phase difference Aip, which when tuned appropriately, 
can flip direction or vanish altogether. Eor comparison, the 
ds - 0 case is included, which has symmetric behavior about 
Aip = 180°. Eor most S layer widths considered, Tx,tot van¬ 
ishes at Aip ss 90°, and Aip =» 270° before reversing direction. 
Since the middle S terminal has ip^ - (P/2, the central S elec¬ 
trode tends to asymmetrically distort the supercurrent about 
Atpn = 180°. Comparing Eig. [^a) with Eig. [^b), it is evi¬ 
dent that the charge supercurrent Jx is not simply correlated 
with the flux of spin current (or equivalently Tx.tot)^ consistent 
with previous worlP. The coupling between ferromagnets is 
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clearly stronger the thinner the S electrodes, where T^^tot is 
larger and tends to change less over a broader range of A(p, 
reflecting a tendency for the magnetization to remain fixed in 
place despite supercurrent variations. Eventually however, for 
sufficient increments in Atf, the net torque will abruptly re¬ 
verse direction. 

The proposed S FS FS system can be considered as a type 
of superconducting magnetic torque transistor, where the flow 
of spin and charge currents are tuned by Atf. This, in turn, 
dictates the torques acting on the exchange fields present in 
the F layers. By minimizing the free energy,® it was shown 
that changes in the supercurrent with respect to relative mag¬ 
netization orientation results in a torque that changes with Atp, 
and vice versa. To underscore the sensitivity of the net torque 
to the phase and relative magnetic orientations. Fig. [^c) il¬ 
lustrates Ti tot as a function Atp for a few orientation angles. 
Pi. When m and h are collinear, i.e., the two exchange field 
alignments in the ferromagnets are parallel {Pi - 0°) or an¬ 
tiparallel {P 2 - 180°) to one another, mxh - 0, and hence 
the net torque is zero (see Eq. ([T0|). The previous P 2 - 90° 
case in (b) is also shown here. For noncollinear magnetiza¬ 
tions, a “static” torque even exists in the absence of a super¬ 
current {Atp - 0). In this case, the effectively inhomogeneous 
magnetization generates a spin current imbalance and torque 
that tends to align the magnetizations. When the magneti¬ 
zations are misaligned, the supercurrent can change both the 
direction and amplitude of the torqud^. In many cases, this 
effect can be attributed to the torque that the equal-spin triplet 
component of the supercurrent (possessing net spin along the 
spin quantization axis) exerts on the magnetization and tends 
to rotate it. As we clearly see from the results presented in 
Fig. I^c), for thin central superconductors with dg < this 
effect can be quite sensitive to the multiple superconducting 
phase differences. 

To investigate further the behavior of the local spin trans¬ 
port, and total torque when varying the ferromagnet orienta¬ 
tion angle Pi, the spatial behavior of spin current through¬ 
out the system is shown in Fig. |^a). The angle Pi de¬ 
scribes the rotation of the in-plane magnetic exchange in Fi: 
hi = /!(0, sinjSi,cosySi). The magnetic exchange field direc¬ 
tion in F 2 does not vary and is directed along z: h 2 - (0,0, h). 
Control of the free-layer magnetization by an external mag¬ 
netic field has experimentally been demonstrated in S IF spin 
valveJ^. The r otation angle can also be manipulated by STT 
switching.SSI We see that in the S region, 5;^ is constant for 
all angles Pi, consistent with the spin-torque conservation law 
[Eq. ([nj] which states that any spatial variations of the spin 
current must generate a torque. The torque thus vanishes in 
the S region, as it should. The F 2 region again exhibits a spa¬ 
tially modulating spin current whose behavior is highly sen¬ 
sitive to the particular orientation angle Pi. The most rapid 
changes in the oscillating S x tends to occur within this fer¬ 
romagnet near the interface with a superconductor. We also 
see that the spin current at the interfaces between the ferro¬ 
magnets and the central S (dashed vertical lines) varies non- 
monotonically, changing sign at Pi - 30°, or vanishing al¬ 
together when the magnetizations are collinear iPi — 0°, or 
Pi — 180°). These observations are consistent with Fig. |^b). 


- 30“ /?j 

- 60" 



xlO^'* 



FIG. 6. (a) The x component of the spin current 5 flowing through¬ 
out a segment of the SfSFS junction, as a function of position 
X. The magnetic exchange orientation in Fi is varied according to 
hi = h(0, sinySi, coSjSi), while in F 2 we have:, hi = {0,0, h). The 
intermediate S width corresponds to ds I^q = \. The legend identi¬ 
fies the different angles Pi used. A current is established via a phase 
difference of Atp = 45° between the outer S layers. Spatial variations 
in are responsible for any torques present in the system. In (b) 
the total torque, Tx.iot, is plotted as a function of Pi. The ds = 0 refer¬ 
ence case is multiplied by a constant factor for comparison purposes. 
Each curve corresponds to a different ds as identified in the legend. 


where the total torque Tx.tot is shown as a function of orien¬ 
tation angle Pi for several S widths. We see that Tx.tot van¬ 
ishes entirely when the two ferromagnets have collinear mag¬ 
netizations, corresponding to parallel {Pi - 0°) or antiparallel 
{Pi - 180°) configurations. The total torque also has the ex¬ 
pected behavior when there is no middle S terminal {ds - 0), 
peaking when Pi x 90°, corresponding to the situation where 
the torque has the greatest tendency to align the magnetic mo¬ 
ments. Including a central S layer is seen to introduce a non¬ 
trivial oscillatory behavior in Ti tot that can cause it to vanish 
(or change direction) multiple times when spanning the full Pi 
range. In effect, the angle Pi that was previously responsible 
for the largest total torque (when ds = 0) is now the angle 
at which there is negligible total torque within the F layers. 
Increasing the S thickness of course reduces the ferromag¬ 
netic coupling and hence reduces the magnitude of the mutual 
torques. The misalignment angle where the maximum torque 
is exerted, ySj, clearly shifts from near the orthogonal configu- 
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central S electrode and into the f 2 layer, generating a spin- 
imbalance. In Fig. |7|a-c), we illustrate the total magnetiza¬ 
tion, tot, in each region as a function of Atp. Here we define 
the quantity OT^. ^ot as the y-component of the magnetization 
spatially integrated over each of the three junction regions of 
interest, and normalized by mo = Using this normal¬ 

ization, the bulk value of the magnetization is equivalent to the 
exchange field value of hi Sf = 0.1. The thicknesses of the Fi, 
F 2 , and S layers are given by dfil^o — 0.1, dril^a - 3.8, and 
ds /^o = 1, respectively. It is evident that for a wide range of 
Aip, a net magnetization exists in each junction region, other 
than when jdx corresponds to relative collinear magnetizations 
{h\ directed along z). In panel (a), my,tot is approximately con¬ 
stant for each j3i, and the overlapping curves at fix - 30°, 150° 
reflect the symmetry about j3x — 90°, where the net magneti¬ 
zation is greatest. Also, the net magnetization in Fx is positive 
for the range of Px shown, since the magnetic exchange inter¬ 
action in the y direction is /ji,y = hsinPx- Due to proximity 
effects, there is an intrinsic net magnetization in the supercon¬ 
ductor that is present even in the absence of current flow. As 
shown in panel (b), the total induced magnetization in the su¬ 
perconductor is finite at Axp - 0° and vanishes at Aip ^ 180°, 
where it switches direction. Finally, in Fig. |^c) the contri¬ 
bution from the magnetization to the total torque observed in 
Fig .|^is evident, where the total magnetization in the wider F 2 
region exhibits the same dependence on Aip as differing 
only in sign. 


IV. CONCLUSIONS 


FIG. 7. The normalized y component of the total magnetization 
as a function of phase, A<p, within the SFSFS system, is 
shown for a few in-plane exchange field orientations, Pi (see legend). 
The Fx region has a thickness corresponding to dfxl^o = 0.1, F 2 has 
dpil^o = 3.8, and the central layer is one coherence length wide, 
i.e., dsl^o = 1. The exchange field has magnitude h = 0.1. Each 
panel (a)-(c) depicts a different region where the net magnetization 
is calculated. The leakage of magnetism into the central 5 region 
[panel (b)] is clearly visible. 


ration (j6* ^ 90°) when ds - 0 towards intermediate magnetic 
configurations corresponding to 60° < < 70°, for ds < ^o- 

To examine the previous behavior of the total STT from a 
different perspective, it is beneficial to recall the simple ex¬ 
pression, Eq. ®, which shows that for a given exchange 
field, the torque arises entirely from the magnetization, m(x). 
Thus, it is insightful to study the details of m(x), which gives 
a measure of the spin polarization in the system responsi¬ 
ble for generating the local spin currents. The out-of-plane 
torque is due to both in-plane components of the magnetiza¬ 
tion: Tx(x) = -(21 fiB)[my(x)h^(x) - hy{x)m^(x)], which clearly 
vanishes outside of the ferromagnet regions where h, - 0. 
The exchange field in the Fx region varies in the yz plane, 
while in F 2 , the only nonzero component to the exchange 
field is h^, so that we have simply, Tx{x) = -{2ffj.B)hz{x)my(x). 
Thus for a mutual torque to exist in the ferromagnets, a y- 
polarized magnetization in Fx must propagate through the 


In conclusion, we have presented a detailed microscopic 
study of the charge and spin supercurrents that can exist in 
SFSFS types of Josephson junction hybrids. The local mag¬ 
netization profiles were then calculated and employed for de¬ 
termining the equilibrium spin transfer torques. We also stud¬ 
ied the associated spin-triplet correlations that arise in these 
hybrids. This was accomplished by solving the BdG equa¬ 
tions over a broad range of geometrical and material parame¬ 
ters. Our investigations revealed how to manipulate and gen¬ 
erate supercurrents with higher order harmonics by varying 
the macroscopic phases in the superconducting electrodes, or 
the relative exchange field orientations. Utilizing the spin con¬ 
servation law, we calculated the spin transfer torque in these 
systems, revealing a number of experimentally viable ways 
in which the magnetization can be controlled in a prescribed 
fashion. Our results demonstrate that, depending on the pa¬ 
rameters considered, these types of ballistic systems can sup¬ 
port supercurrents that can be tuned to contain primarily the 
first or second harmonics in the current-phase relations. We 
studied the 7r-27T harmonic crossovers and determined the ex¬ 
perimentally desirable conditions in which to reveal the sec¬ 
ond harmonic supercurrents in these systems. We also showed 
that the equilibrium spin transfer torques can be well con¬ 
trolled by simply modulating the macroscopic phases of the 
three S electrodes in addition to the other system parameters 
such as the sizes of the ferromagnets and central supercon¬ 
ductor electrodes, or the relative magnetization alignments. 
















9 


These findings are suggestive of a phase-tunable supercon¬ 
ducting transistor based on STT switching. 
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Appendix A: Numerical procedure for solving the BdG 
equations 


The numerical procedure used in calculating the spin and 
charge currents involves first expandinj^ the quasiparticle 
amplitudes in terms of a complete set of N basis functions; 

M sm(kgX)if/g(kg), (A 1 ) 

g=0 

where we define tf/nix) = (m„i-(x), y„i-(v), v„|(x)), and 

^?i)' The wavevector kq - qnjd is dis¬ 
cretized in terms of the system width d, taken to be large 
enough so that the results become independent of d. The next 
step involves Fourier transforming the real-space BdG equa¬ 
tions (Eq. Q), resulting in the following set of coupled equa¬ 
tions in momentum space: 


1 

■ 

1 

• 

+ 

0 

A 

'iij 


'uj 

1 

1 

Hq + 

A 

0 

Mi 


Mi 

0 

A* 

1 

1 

-hx - ihy 



V>| 

, A* 

0 

—hx + ihy 

-(Ho + K)^ 

Ti. 


Ti, 


(A2) 


Here we have defined Mq- = (mio-, M 20 -, • ■ ■, mvo-), Vo- = 
V 2 cr, ■ ■ ■, y^o-), and the matrix elements, 


;S— + f ± - F 
Zm 


2 r 

Ho(q,q) = j dx 
d Jo 

2 

A(q,q)-- I dxA{x)^\\\{kqX)m\{kq,x), 

d Jo 


s,m{kqx)?,m{kqfx), (A3) 


(A4) 


hiiq 




dx hiix) ^\n{kqx) ^\n{kq'X), i — x, y, z. 


(A5) 


Our numerical procedure for calculating the supercurrent in¬ 
volves assuming a constant amplitude and phase for the pair 
potential in each S layer, thus providing the physically neces¬ 
sary source or sink of current, via the external electrodes. We 
then expand the pair potential via Eq. ( |A4| l. Similarly the ex¬ 
change field and free particle Hamiltonian are expanded using 
Eq. ( |A5| l and Eq. ( |A3| l respectively. We then find the quasipar¬ 
ticle energies and amplitudes by diagonalizing the resultant 
momentum-space matrix (Eq. Once the momentum- 

space wavefunctions and energies are found, they are trans¬ 
formed back into real-space via Eq. (All, and the currents 
and magnetic moments are calculated as described in Sec.[n| 
As mentioned above, the current source arises from the non 
self-consistent region where we take A(.r) to be a piecewise 


constant with prescribed macroscopic phases in the S elec¬ 
trodes. The widths of the two outer S terminals are sufficiently 
large {d^ » ^0) so that the system boundaries have a negligi¬ 
ble influence on the results. By taking the divergence of the 
current in Eq. (j^ and using the BdG equations (Eq. ([^), we 
find, 

= 2eIm|A(x)^ [ul^v„i + tanh(^) 

V n 

(A6) 

where the terms in the summation constitute the usual self- 
consistency equatiorPS for A(x). Thus, only when self- 
consistency in A(x) is achieved, does the right hand side 
of Eq. ( |A6| l vanish, and current is conserved. If the self- 
consistency condition is not strictly satisfied, the terms on the 
right act effectively as sources of current, except of course 
within the ferromagnet regions, where A(x) = 0. 


Appendix B: Spin current and spin transfer torque 

The spin current can be found by using the Heisenberg pic¬ 
ture. Eirst we determine the time evolution of the spin density, 
V(x), 


-{r](x))^i{m,v(x)]), (Bl) 

where the effective BCS Hamiltonian 'Tf is written. 


•7/ = 



(x)['77o(x) - h(x) ■ cr]i/r(x) 


-I- A(x)i/fj(x)(/r|(x) -I- A*ix)ifri{x)i//^{x)]. 


(B2) 


Here we define, i/^(x) = {ij/Jx), i/f|(x))^, and the spin density 
operator J7(x): ri{x) - (x)cr\l/{x). Inserting the Hamiltonian 

(Eq. (IB^), into Eq. ®, we end up with the spin continuity 
equation: 


— {'nix)) + -^Six) = t{x), (B3) 

at ox 

where the spin transfer torque t is written in terms of the ex¬ 
pectation value, t{x) - 2{\l/\x){cry.h\xli{x)). Using the fact 
that the spin density is simply related to the magnetization m 
via m(x) - -/jb inix)), we end up with Eq. ([Tgi. Similarly, 
the spin current S is given by: 

S(x) = -^|iA'(x)cr^^^/r(x)j - ^^t/r^(x)jcri/r(x)|. (B4) 

Lastly, we insert the Bogolibuov transformations, (/'t(x) = 
Ijn(unJx)y„-vl^(x)yl), and i/^x(x) = 2„(M„|(x)y„ -r v*^(x)yl), 
and use conventional rules for the thermal averages: (ylym) - 
Snmfn, {jmyl) = 5„„(1 - /„), and {ynjm) = 0, to arrive 
at Eqs. 0-0- Note that 5o-(Ji:) represents the spin current 
flow along the x direction in configuration-space, with indices 
cr = X, y, z in spin-space. 
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